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ABSTRACT 

We  provide  a  map  0  tq  which  associates  each  finite  subset  0  C  C*  with  a  polynomial  space 
TT©  from  which  interpolation  to  arbitrary  data  given  at  the  points  in  0  is  possible  and  uniquely 
so.  Among  all  polynomial  spaces  Q  from  which  interpolation  at  0  is  uniquely  possible,  our  ;r© 
is  of  smallest  degree.  It  is  also  D-  and  scale-invariant.  Our  map  is  monotone,  thus  providing  a 
Newton  form  for  the  resulting  interpolant.  Our  map  is  also  continuous  within  reason,  allowing  us  to 
interpret  certain  cases  of  coalescence  as  Hermite  interpolation.  In  fact,  our  map  can  be  extended  to 
the  case  where,  with  each  d  €  0,  there  is  associated  a  polynomial  space  Pg,  and,  for  given  smooth 
/,  a  polynomial  q  &  Q  is  sought  for  which 

p(D)(/-q)(0)  =  O,  VpePfi,  ^6  0. 

We  obtain  itq  ais  the  “scaled  limit  at  the  origin”  (exp©);  of  the  exponential  space  exp©  with 
frequencies  0,  and  base  our  results  on  a  study  of  the  map  H  defined  on  subspaces  H  of 

the  space  of  functions  analytic  at  the  origin.  This  study  also  allows  us  to  determine  the  local 
approximation  order  from  such  H  and  provides  an  algorithm  for  the  construction  oi  from  any 
basis  for  H. 


AMS  (MOS)  Subject  Classifications:  primary  41A05,  41A63,  41A10;  secondary  65D05,  41A30 
Key  Words:  exponentials,  polynomials,  multivariate,  interpolation,  Newton  form 


^  supported  in  part  by  the  National  Science  Foundation  under  Grant  No.  DMS-8701275  and  by 
the  United  Slates  Army  under  Contract  No.  DAAL03-87-K-0030 


UNIVERSITY  OF  WISCONSIN-MADISON 
CENTER  FOR  THE  MATHEMATICAL  SCIENCES 


On  multivariate  polynomial  interpolation 
Carl  de  Boor^  k  Amos  Ron 


ABSTRACT 

We  provide  a  map  0  jt©  which  associates  each  finite  subset  0  C  C*  with  a  polynomial  space 
TT©  from  which  interpolation  to  arbitrary  data  given  at  the  points  in  0  is  possible  and  uniquely 
so.  Among  all  polynomial  spaces  Q  from  which  interpolation  at  0  is  uniquely  possible,  our  tt© 
is  of  smallest  degree.  It  is  also  D-  and  scale-invariant.  Our  map  is  monotone,  thus  providing  a 
Newton  form  for  the  resulting  interpolant.  Our  map  is  also  continuous  within  reason,  allowing  us  to 
interpret  certain  cases  of  coalescence  as  Hermite  interpolation.  In  fact,  our  map  can  be  extended  to 
the  case  where,  with  each  ^6  0,  there  is  associated  a  polynomial  space  Pg,  and,  for  given  smooth 
/,  a  polynomial  g  6  Q  is  sought  for  which 

p{D)(f  -  q)ie)  =  0,  Vp€P«,  <?€©. 

We  obtain  ttq  as  the  “scaled  limit  at  the  origin”  (exp©)!  of  the  exponential  space  exp©  with 
frequencies  0,  and  base  our  results  on  a  study  of  the  map  H  ^  defined  on  subspaces  H  of 
the  space  of  functions  analytic  at  the  origin.  This  study  also  allows  us  to  determine  the  local 
approximation  order  from  such  H  and  provides  an  algorithm  for  the  construction  of  Hi  from  any 
basis  for  H. 


AMS  (MOS)  Subject  Classifications:  primary  41A05,  41A63,  41A10;  secondary  65D05,  41A30 
Key  Words:  exponentials,  polynomials,  multivariate,  interpolation,  Newton  form 


Authors’  affiliation  and  address: 
Center  for  Mathematical  Sciences 
University  of  Wisconsin-Madison 
610  Walnut  St. 

Madison  WI  53705 


^  supported  in  part  by  the  National  Science  Foundation  under  Grant  No.  DMS-8701275  and  by 
the  United  States  Army  under  Contract  No.  DAAL03-87-K  0030 


On  multivariate  polynomial  interpolation 
Ccirl  de  Boor  &  Amos  Ron 

r. 

1.  Introduction  /  ^ 

/  \  ■ 

\  / 

The  gener^li«S[tion  of  upivariate  polynomial  interpolation  to  the  multivariate  context  is  made 

difficult  byAhe  fact  that  one  has  to  decide  just  which  of  the  many  of  its  pice  properties  to  preserve, 
as  it  iadmpossible  to  preserve  them  aU.  Particirlarly  annoying  is  the  fact  that  the  dimeimions  of 
stand^d  multivariate  polynomial  spaces,  such  as  tt*,  make  up  only, ai  small  subset  of hence 
w^annot  hope  to  interpolate  uniquely  at  an  arbitrary  pointset  •  &  (t  C“  from  an  appropriate 

Further,  even  when  we  have  dimTr^  points  at  hand,  they  may  fail  to  be  total  for  lei'^ence 
interpolation  at  these  points  from  ^  may  still  not  be  possible,  (let )  ^ - 

For  these  reasons,  generalizations  have  stressed  some  aspects  of  polynomial  interpolation  and 
ignored  others.  For  example,  there  are  various  efforts  (see,  e.g.,  [CY],  [GM])  to  identify  finite 
sets  0  for  which  it  is  easy  to  construct  polynomial  Lagrange  functions,  i.e.,  polynomials  pe  with 
pe(T)  =  6e,T-  Except  for  special  circumstances,  it  is  usually  hard  to  ascertain  the  degree  of  the 
resulting  interpolant  or  the  maximal  k  for  which  iCk  lies  in  the  range  of  this  interpolant.  A  totally 
different  effort,  associated  with  the  name  Kergin  (see  [K],  [M]),  retains  the  fact  that,  with  an 
arbitrary  set  0  of  cardinality  A:  +  1,  we  interpolate  from  Xfc.  The  additional  degrees  of  freedom 
available  in  a  multivariate  context  Kergin  disposes  of  in  such  a  way  that  ‘natural’  meanvalue 
theorems  continue  to  hold. 

In  this  paper,  we  take  a  different  tack.  Given  any  finite  set  0,  we  determine  a  corresponding 
polynomial  space  ttq  from  which  interpolation  at  0  is  ‘correct’,  i.e.,  is  possible  and  uniquely  so. 

We  show  that  itq  is  translation-  and  scale-invariant,  and  that  it  is  a  polynomial  space  of  least 
degree  from  which  interpolation  at  0  is  correct.  We  also  show  that  the  resulting  map  0  i-»  x© 
is  monotone  (as  a  map  from  sets  to  sets),  making  it  natural  to  introduce  a  Newton  form  for  the 
resulting  interpolant.  Further,  we  show  that  the  map  can  be  extended  in  a  natural  way  to  Hermite 
interpolation,  where  we  allow  some  of  the  B  to  coalesce. 

In  fact,  given  arbitrary  finite-dimensional  polynomial  spaces  we  provide  such  a  polynomial 
space  of  least  degree  from  which  “generalized  Birkhoff-Hermite”  interpolation  is  correct,  i.e.,  over 
which  the  linear  space  spanned  by  the  linear  functionals  of  the  form  [B]p{D),  with  p  G  Pg  and 
9  £  Q,  is  minimally  total.  Here,  [9]f  ;=  f{9),  and,  to  recall,  a  space  A  of  linear  functionals  is  total 
for  H  if  the  only  h  £  H  foi  which  Ah  =  0  for  all  A  €  A  is  /i  =  0. 

The  following  notation  and  terminology  will  be  used  throughout.  The  collection  of  all  poly¬ 
nomials  on  C‘  (or  whatever  other  space  the  context  might  indicate)  is  denoted  by  tt;  TTk  denotes 
the  collection  of  all  those  polynomials  of  (total)  degree  <  k,  i.e.,  Wk  :=  span  (O")!^!^^,  with 
0“  ;  I  (-<■  We  also  find  it  convenient  to  use  7r<jt  for  the  space  of  polynomials  of  degree  <  k. 

For  any  p  6  x,  we  denote  by  p{D)  the  corresponding  constant  coefficient  differential  operator;  in 
particular,  D°‘  :=  differentiation  with  respect  to  the  jth  argument.  We  make 

good  use  of  the  representation  of  the  linear  functional  {9]p{D)  on  tt  as  9  p*E^q  =  q*{egp),  with 

q-p  :=  (g(Z))p)(0)  =  j;Z?“?(0)i?>(0)/a!, 

a 

with  eg  :  X  >-*  exp{{9,x)),  and  with  E  tlie  shift,  i.«*,  F*/  ■=  /(•  4-  B).  We  also  use  the  “least  term” 
of  a  function  /  analytic  at  the  origin,  i.e.,  the  homogeneous  polynomial  fi  of  largest  degree  j  for 
which  /(z)  =  /i(z)  -t-  o()|ilp)  as  z  — ►  0.  This  notion  makes  (formal)  sense  on  the  larger  space  of 
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formal  power  series,  as  does  the  notion  of  a  Z)-invariant  subspace,  i.e.,  a  subspace  invariant  under 
differentiation. 

The  paper  is  laid  out  as  follows.  After  a  short  discussion  of  basic  properties  of  the  “leading 
term”  and  the  “least  term”  fi  of  /,  we  show  in  Section  3  that,  for  any  finite- dimensional  space 
H  of  sufficiently  smooth  functions,  Hy  :=  span{/x  :  /  G  is  a  scale-invariant  polynomial  space 
of  the  same  dimension  as  H  which  determines  the  “local  approximation  order”  from  H.  We  also 
establish  other  properties  of  Hi,  such  as  the  fact  that  p‘\(D)Hi  C  {p{D)E)i,  of  use  later,  and 
explore  Hi  for  the  space  H  =  exp0  spanned  by  the  exponentials  eg,  9  £  0.  Section  4  is  devoted  to 
the  important  observation  that,  among  aU  polynomial  spaces  P  for  which  the  corresponding  space 
P  :=  :  p  G  P}  of  linear  functionals  is  minimally  total  for  H,  Hi  is  of  least  degree  in  the 

sense  that,  for  all  j,  dim(P  Cl  Xj)  <  dim(P|  fl  Xj).  The  resulting  linear  projector  given  by  H  and 
11*1  is  exploited  in  Section  5  in  the  derivation  of  an  algorithm  for  the  construction  of  Hi  from  any 
basis  for  H.  The  dependence  of  Hi  on  H  is  explored  in  Section  6;  the  main  result  is  that  the  map 
H  Hi  is  continuous  at  H  if  and  only  if,  for  some  m,  x<„i  C  c  x„,,  a  property  of  H  which 
we  term  “regular”,  for  want  of  a  better  word. 

The  remainder  of  the  paper  is  devoted  to  the  specific  choice  H  =  of  exponentials. 

The  fact  that  its  least  part  Hi  supplies  correct  conditions  for  interpolation  from  H  is  used  in 
Section  7  to  conclude  by  duality,  as  in  [DR],  that  it  is  possible  to  interpolate,  and  uniquely  so,  from 
the  polynomial  space  Hi  using  the  interpolation  conditions  [6]p{D),p  G  Pe,9  G  0.  The  special 
case  Pg  =  xo,  aU  9,  leads  in  Section  8  to  Lagrange  interpolation  from  x©  :=  (exp©)x,  witn  the 
algorithm  from  Section  5  providing  information  needed  for  the  Newton  form  for  the  interpolant. 
The  connection  between  coalescence  of  such  interpolation  points  and  osculatory  interpolation  is 
explored  in  the  final  section. 

In  a  subsequent  paper,  we  verify  that  various  forms  of  multivariate  Lagrange  interpolation  now 
in  the  literature  are  special  cases  of  the  scheme  proposed  here.  In  a  diflFerent  paper,  we  use  Hi  to 
simplify  and  extend  results  from  box  spline  and  exponential  box  spline  theory. 

2.  The  least  term  of  an  analytic  function 

We  denote  by  p^  the  leading  term  of  the  polynomial  p.  For  p  ^  0,  this  is  the  (unique) 
homogeneous  polynomial  for  which 


deg(p-pt)  <  degp. 

For  completeness,  we  take  the  zero  polynomial  to  be  its  own  leading  term.  We  note  'that  {pq)i  = 
Pt9t- 

We  also  single  out  the  least  term  /;  (read  ‘/  least’)  of  a  polynomial  or,  more  generally,  a 
function  /  analytic  at  the  origin,  and  mean  by  this 

/i  :=  Tif, 

with  j  the  smallest  integer  for  which  Tjf  0,  with  Tj/  the  Taylor  polynomial  of  degree  <  j  for  / 
at  the  origin,  i.e., 

(2.1)  T^/  :=  Y.  iTD°m. 

I“l<i 

and  with  0“  :  i  i->  x“/a!  the  normalized  power  function.  For  completeness,  we  take  the  zero 
function  to  be  its  own  least  term.  We  note  that 

(2.2)  ifg)i  =  fi9i 
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and  that,  for  any  invertible  matrix  A, 


(2.3) 


{foA)i  =  fioA. 


We  are  interested  in  fi  because  it  describes  the  behavior  of  /  near  the  origin.  Precisely,  /|  is 
the  homogeneous  polynomial  of  largest  degree  j  for  which 

(2.4)  /(*)  = /i(a:)  +  o(||i|p)  as  i -h- 0. 

Consequently,  with  j  :=  deg/^, 

(2.5)  /i  =  lim/(t-)/f^ 


in  the  pointwise  sense,  as  follows  readily  from  L’Hopital’s  rule.  In  fact,  (2.5)  remains  true  if  we 
take  the  limit  in  the  sense  of  formal  power  series,  i.e.,  in  the  sense  that,  for  all  a, 

[0]Z?“/i  =  limj0](J9“/)(i-)A^ 


3.  The  limit  at  the  origin  of  a  space  of  functions  analytic  there 
In  this  section,  we  consider  subspaces  of  the  space 

Ao 

of  all  functions  analytic  at  the  origin,  with  the  topology  of  formal  power  series.  For  any  subspace 
H  of  i4o,  we  consider  its  “limit  at  the  origin”,  i.e.  (with  (2.5)),  the  polynomial  space 

(3.1)  Hi  :=  span{/|  :  /  6  H}. 

We  note  that  Hi  is  scale-invariant  since  it  is  spanned  by  homogeneous  polynomials. 

We  were  led  to  in  the  analysis  of  the  local  approximation  order  from  H.  By  definition, 
this  is  the  largest  integer  d  for  which,  for  every  /  G  C'°°(ni*),  there  exists  h  £  H  so  that 

(/-  /i)(x)  =  0(lli!|’^)  as  I  0. 

The  following  lemma  is  of  use  in  the  discussion  of  approximation  order. 

(3.2)  Lemma.  If  7r<fc  C  E^,  then  there  exists  a  continuous  linear  projector  Tk,H  on  Aq  into  E 
for  which  TkTk,H  =  Tk. 

Proof:  If  7r</t  C  E i,  then,  for  each  |q|  <  k,  there  exists  fa  £  E  with  {fa)i  =  follows 

that  the  matrix  (f^“//3(0))|o,|,|/3|<fc  is  unit  triangular,  hence  invertible.  This  implies  that  we  can 

find  (fifa)|Q|<fc  in  E  dual  to  (/  ^"/(0))|a|<fc  (i-6->  satisfying  D°‘gp{0)  =  dap),  and  this  implies 

that,  for  each  /  £  Aq,  h  :=  Tk,^/  ■=  23affa^“/(0)  i®  ^  satisfies  Tkh  ~  Tkf.  ♦ 
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(3.3)  Corollary.  The  local  approximation  order  of  a  finite-dimensional  subspace  H  of  Ao  equals 
the  largest  integer  d  for  which  7r<d  C  JTi- 

Proof:  Let  d  be  the  local  approximation  order  from 

Having  {f  —  h){x)  =  0(11x11*^)  as  i  -+  0  is  the  same  as  having  deg{f  -h)i  >  d.  If,  in  particulzir, 
/  G  then  this  cam  only  happen  if  =  fi-  Since  7r<d  =  (7r<(f)^,  this  shows  that  7r<d  c  Hi- 
Conversely,  if  3r<fc  C  Hi,  then,  by  the  lemma,  there  is  a  linear  map  Tk,H  into  H  with  7^(1  - 
Tk,H)  =  0.  This  implies  that,  for  any  /  6  Ao,  h  :=  Tk,Hf  is  in  H  and  satisfies  (/-/i)(x)  =  0(||x||*'), 
hence  d>  k.  4 

Further  study  of  Hi  led  us  to  the  results  on  polynomial  interpolation  to  be  detailed  in  subse¬ 
quent  sections.  In  preparation,  we  now  discuss  various  properties  of  Hi- 

Let  H  be  a  finite-dimensional  subspace  of  Ao.  We  observe  that,  for  f  €  H,  deg  fi  =  j  if  and 
only  if  /  6  (ker^rTj)\(kerH  T,+x),i.e.,  if  and  only  if  fi  ^  7)^2(ker//7j‘)\0,  with 

(3.4)  keiff  Tj  :=  ker(Tj|H)- 
Since  also 

H  =  ker/f  To  D  ker//  Ti  D  kerjf  Tj  5  •  •  •  > 


hence 

dim  Tj+i  (kerjf  Tj )  =  dim  kerif  Tj  —  dim  ker^ 


we  conclude  that 


oo  oo 

Hi  =  ^T,+i(kerHT,)  =  0T,+i(kerHr,) 
i=o  i=o 


and 


k 

dim(H;  n  :rfc)  =  dim ker/f  Tj  —  dimker// Tj+i  =  dimker/fTo  —  dimker^Tfe+i  =  dimrfc+i(i5'). 

j=0 


We  have  proved; 

(3.5)  Proposition.  Hi  is  a  scale-invariant  space  of  polynomials  of  the  same  dimension  as  H.  In 
fact,  for  every  j, 

(3.6)  dim(5'i  n  x^j)  =  dimTj{H). 

Also,  (Hi)i  =  Hi,  and  (TjH)i  =  Hi  for  all  sufficiently  large  j,  and  H  =  Hi  in  case  H  is  a 
scale-invariant  polynomial  space. 

For  the  particular  space  H  ;=  span  (l  -|-  ()^’®,()”’^),  one  computes  that  Hi  =  span  (l,()°4)  ^ 
H  =  T2{H),  thus  illustrating  that  Tj{H)  and  D  need  not  be  equal  in  general  (even  though 
they  are  always  of  the  same  dimension). 

Next,  consider  the  effect  of  multiplying  all  the  elements  of  H  by  some  /  G  Ao,  i.e.,  the 
relationship  between  H  and 

fH  :=  {fg\  g  G  H). 

We  deduce  from  (2.2)  the  following  observation. 


4 


(3.7)  Proposition.  For  any  f  ^  Aq  sa,tisfying  /(O)  0, 

{fH)i  =  Hi. 

The  interaction  of  diflferentiation  with  the  map  ^  is  determined  by  the  fact  that,  for 

any  p  £  tt  and  any  /  G  Aq, 

(3.8)  p{D)f  =  p-^{D)fi  +  terms  of  higher  degree. 

This  implies  that  {p{D)f)i  =  pi(D)fi  in  case  p^{D)fi  0  and  so  proves  the  following. 

(3.9)  Proposition.  For  every  p 

(3.10)  C  {p{D)H)^. 

(3.11)  Corollary.  Ifp{D)  annihilates  H ,  then  Pj(.D)  annihilates  Hi- 

(3.12)  Corollary.  If  H  is  D-invariant,  then  so  is  Hi- 

Proof:  For  every  y  6  IR*,  we  have  Dy{Hi)  C  (DyH)i  by  (3.9)Proposition,  while  (DyH)i  C 
since  DyH  C  H  hy  assumption.  4 

If  H  consists  of  functions  analytic  on  some  domain  <?,  then  it  makes  sense  to  consider  the 
“limit  of  H  at  z"  for  any  z  ^  G.  If  /f  is  Zl-invariant,  hence  translation-invariant,  we  expect  aU 
these  limits  to  coincide.  The  Corollary  confirms  this  expectation. 

As  an  example,  consider  the  space 

(3.13)  exp©  :=  span{ee  :  0  €  0} 
of  simple  e.xponentials  with  frequencies  0.  Here  and  below, 

e«  :  1 1-*  exp((0,x}). 

Since  exp©  is  Z?-invariant,  its  limit  at  any  point  is  just  (exp©)|.  For  its  construction,  we  can 
proceed  as  follows.  Define 

PaJ  :=  ^(^•)^«(^). 

«e© 

with 

(0-y  :  X  >->■  {9,xy  fjl. 

Then 

(3.14)  TT©  :=  (exp©)!  =  0{pa.j  :  Pa,i  =  0  for  i  <  j}. 

i 

The  D-invariance  of  tt©  can  also  be  seen  directly,  as  follows:  For  y  G  IR*,  DyPaj  =  Payj-i,  with 
ay  :=  {a{d){9,y}) Hence  ifpo,,-  =  0  for  all  i  <  j,  then  Pay,i-i  =  DyPa,i  =  0  for  all  i  -  1  <  j  -  1, 
therefore  DyPaj  G  Jr©  if  paj  G  ^0. 

4.  The  dual  of  H 

As  it  turns  out,  the  construction  of  can  be  carried  out  by  a  bootstrap  procedure  which 
uses  interpolation  from  H.  For  this  reason  (and  others),  we  now  show  that  the  dual  of  H  can  be 
represented  by  Hi- 

Abstractly,  interpolation  from  H  can  be  described  as  the  task  of  finding,  for  given  /  G  Ao,  an 
h  ^  H  fox  which  Xh  =  Xf  for  all  A  in  some  linear  space  A  of  linear  functionals  on  Aq.  VVe  call  A  the 
(space  of)  interpolation  conditions  for  this  particular  interpolation  problem.  We  call  the  problem 
correct  if  there  is,  for  each  /,  exactly  one  solution  h. 

For  completeness,  we  recall  (without  proof)  the  following  well  known  characterizations  of 
correctness. 
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(4.1)  Lemma.  Let  H  and  A  be  Haite-dimensional  linear  subspaces  of  a  linear  space  X  (overC) 
and  Its  dual,  respectively.  Then  the  following  are  equivalent. 

(I)  The  Interpolation  problem  given  by  H  and  A  Is  correct. 

(II)  With  (Aj)i  any  basis  for  A,  the  Unear  map  E  —>■  :  h  i~*  (^jh)^  Is  cne-one  and  onto. 

(III)  A  Is  minimally  total  for  H. 

(Iv)  A  can  be  used  to  represent  the  dual  H'  of  H  In  the  sense  that  the  map  A  -*  E'  :  A  A|//  Is 
one-one  and  onto. 

If  the  interpolation  problem  given  by  E  and  A  is  correct,  then  it  defines  a  linear  projector, 
P  :=  Ph,a  say,  by  the  rule  that,  for  any  f,Pf^E  and  AP/  =  A/  for  all  A  6  A.  We  will  be 
interested  later  in  the  dependence  of  Ph,k  OJI  siod  A.  For  this,  we  remark  that  P  can  be  written 
in  the  form 

P  =  V{MV)-'^M, 

with  V  any  coordinate  map  for  E  (i.e.,  V  :C"*  -*  E  '.  a  ^  for  some  basis  (t>j)  for  E), 

and  M  the  dual  of  any  coordinate  map  for  A  (i.e.,  M  •.  X  —*  C"  :  /  (Aj/)  for  some  basis  (\j) 

for  A).  This  implies  that  P^.a  is  close  to  Ph'A'  oase  E'  and  A'  have  bases  close  to  bases  of  E 
and  A  respectively. 

Concretely,  we  are  interested  in  using  linear  functionals  of  the  form 

(4.2)  r  :  /  -  {piD)f)iO)  =  ^:FMJ5“/(0)/a! 

Ct 

with  p  G  IT.  These  are  continuous  linear  functionals  on  Aq  and  even  on  C*'(0)  (over  the  complex 
scalars).  The  map  pt-^p*  is  skew-linear  and  one-one,  hence  provides  a  skew-linear  embedding  of  tt 
in  the  dual  of  Aq. 

(4.3)  Proposition.  For  any  Unlte-dimenslonal  linear  subspace  E  of  Aq,  the  linear  space  Is 
minimally  total  for  E. 

Proof:  For  any  /  €  E\0,  p  :=  /j,  G  Ei  and  p*f  =  p*p  >  0.  This  implies  that  the  only 
f  €  E  with  p*/  =  0  for  all  p  €  is  /  =  0,  i.e.,  E^  is  total  for  E.  On  the  other  hand,  since 
dim  E  ^  =  dim  E^  =  dim  E,  no  proper  subspace  of  E^  could  be  total  for  E.  4k 

We  conclude  that  E  ^  can  be  used  to  represent  the  dual  of  E. 

Of  course,  there  exist  polynomial  spaces  other  than  E^  that  provide  correct  interpolation 
conditions  for  E.  The  next  result  shows  that,  compared  to  all  these  spaces  P,  E^  is  of  least 
degree  in  the  sense  that 

(4.4)  dim(P'i^  D  itj)  >  dim(P  n  tTj)  Vj. 

(4.5)  Theorem.  Among  all  polynomial  spaces  P  for  which  T'  Is  minimally  total  for  E,  E^  Is  of 
least  degree. 

Proof:  Let  P  be  a  basis  for  P  D  itj.  Since  P*  is  minimally  total  for  E,  the  sequence  'S'  must 
be  linearly  independent  over  E .  On  the  other  hand,  for  amy  p  £  'Kj,p‘  =  p*Tj+i.  Hence,  B  must 
already  be  linearly  independent  over  Tj+i{E).  Therefore,  with  (3.6), 

dim(Pi  n  TTj)  =  dimTj+i(E)  >  #P  =  dim(P  D  Tj). 


4k 
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5.  The  construction  of  Hi 


Denote  by  Tji  the  linear  projector  given  by  H  and  H  i.e.,  the  linear  map  characterized  by  the 
fact  that  its  range  is  H  and  that  p^Taf  =  p*  f  fot  all  p  £  H i-  The  following  strong  monotonicity 
property  will  be  of  use  in  the  construction  of  Hi  for  specific  choices  of  H. 

(5.1)  Proposition.  If  K  =  H  +  span{/i},  then  Ki  =  Hi  +  span{(/i  -  Tf{h)i}. 

Proof;  There  is  nothing  to  prove  in  case  h  e  H,so  assume  that  h  ^  H,  hence  k  :=  h-  Tnh  ^  0, 
therefore  also  ft  0,  while  £  Ki.  Since  Hi  C  Ki  and  dim/f^  =  dim/f  =  dim^  +  1  = 
dim +  1,  we  therefore  need  only  show  that  ki  ^  Hi. 

For  this,  observe  that,  by  construction,  ^A:  =  0  for  all  p  G  Hi.  Hence,  from  (3.8),  fki  =  0  for 
all  homogeneous  p  G  Hi.  If  now  ki  G  Hi,  then,  in  particular,  kiki  =  0,  i.e.,  Ar^  =  0,  a  contradiction. 

4 

This  proposition  suggests  a  simple  Gram-Schmidt-like  algorithm  for  the  conversion  of  a  basis 
(pj)  for  H  into  a  basis  (<7^)  for  H  to  which  (r^)  :=  (gj^)  is  bi-orthogonal  with  respect  to  the 
(complex)  pairing 

(,)  :  IT  X  Ao  -*  C  :ip,q)>-*  p*q, 

hence  provides  the  homogeneous  orthogonal  basis  (rj)  for  Hi.  The  idea  is  simple.  Assume  that  we 
have  already  determined  such  a  basis  {qj)j<k  for  Hk  :=  span  {pj)  Then  we  compute 

qk  :=  (1  -  TWJPfe  =  Pfc  - 

and  note  that  qk  0  since  pk  ^  Hk  -  ran  T/f* .  Consequently, 

rfc  :=  qki 

is  not  zero,  and  {Tk,qk)  0.  By  (5.1)Proposition,  we  know  that  (.fffc+i)i  =  span  Also,  by 

construction, 

(5.2)  {rj,qk)  =  0  j<k. 

But  there  is,  off-hand,  no  reason  to  expect  that  {rk,qj}  =  0  for  j  <  k.  Of  course,  if  deg  r*  <  deg  rj, 
then  we  have  {rk,qj)  =  0  trivially.  If  degr*  =  deg  rj,  then  (rk,qj)  =  {rj,qk)  =  0.  But  for 
degrfc  >  degfj,  we  may  well  have  {rk,qj)  ft:  0.  In  that  case,  we  simply  modify  qj  appropriately, 
setting 

(5.3)  qj  :=  qj  -  qk  if  deg  r*  >  deg  . 

By  (5.2),  this  does  not  change  the  bi-orthogonality  of  (rj)j<jt  and  {qj)j<k-  Also,  since  it  modifies 
qj  at  terms  of  order  degr*  and  higher,  it  does  not  change  qj^,  i.e.,  it  does  not  change  the  fact  that 
rj  =  for  j  <  k.  In  this  way,  we  have  now  at  hand  a  basis  of  the  promised  sort  for  Hk+i- 
For  easy  reference,  we  collect  the  result  of  the  last  paragraph  in  the  following. 


(5.4)  Algorithm.  Given  the  basis  [pj)  of  the  finite-dimensional  subspace  H  of  Aq. 

For  k  =  1,2,...,  carry  out  the  following  three  steps: 

Step  2.  r*  qki 

Step  3.  qj  *-  qj  -  qk  if  degr*  >  degr^-. 

{rk,qk) 

Then  (rj)  is  bi-orthogonal  to  (qjJ  and  provides  a  homogeneous  orthogonal  basis  for  Hi. 

For  the  calculations,  it  is  useful  to  observe  that  only  inner  products  with  the  homogeneous 
polynomials  rj  are  required.  This  means,  in  particular,  that  the  calculation  could  be  carried  out 
with  TmPk  rather  than  pk,  for  some  m  which  is  determinable  a  priori  in  case  H  is  D-invariant. 

Numerically,  the  calculation  is  challenging  only  because  it  requires  the  determination  of  the 
least  part  rk  of  qk-  When  using  finite-precision  arithmetic,  it  may  be  necessary  to  replace  ‘least 
part’  by  ‘significant  least  part’  in  order  to  avoid  use  of  a  least  part  that  turned  out  not  to  be  zero 
only  because  of  the  noise  in  the  calculation.  Concretely,  this  means  that  one  takes  rjt  to  be  the 
homogeneous  part  of  qk  of  lowest  degree  which  is  not  significantly  smaller  than  the  corresponding 
part  of  Pk.  Considerations  of  this  kind  could  be  used  to  establish  under  what  circumstances  Hi 
depends  continuously  on  H .  In  the  next  section,  we  choose  to  settle  this  question  by  a  more  direct 
route. 


6.  Continuity  of  the  map  H  Hi 

In  the  discussion  later  of  Hermite  interpolation  as  the  limit  of  Lagrange  interpolation,  it  will 
be  important  to  understand  to  what  an  extent  Hi  depends  continuously  on  H.  Precisely,  we  wish 
to  know  under  what  circumstances  limt_o  =  {limf_»o  H{t)) Since  Aq  is  a  metric  linear 

space,  we  use  the  gap  between  subspaces  as  a  means  of  defining  the  statement  H  =  limt— o  H{t). 
For  this,  recall  the  standard  definition  of  the  gap 

gap{H,K)  :=  max{dist(fl^  D  5,/<'),dist(/$'  r\B,H)} 

between  subspaces  H,K  of  the  metric  linear  space  X.  Here,  dist(y,  Z)  :=  supj,gy-inf*g^dist(y,z), 
and  B  :=  Hi(0)  :=  {i  G  A"  :  dist(i,0)  <  1}. 

If  H  is,  in  particular,  a  finite-dimensional  subspace,  then  H  =  limH^(t)  iff,  for  some  (every) 
basis  {hj)  of  H  and  all  small  enough  t,  there  is  a  corresponding  ba.sis  {hj{t))  for  H{t)  for  which 
hj  =  limt_o  hjit)  for  all  j.  For  example,  by  (2.5), 

(6.1)  lini{/i(t-)  -.h^H)  =  Hi. 


As  an  illustration  of  the  possible  lack  of  continuity  in  the  map  H  Hi,  consider  H(t)  =  expQ(j) 
with  s  =  2  and  Q{t)  :=  {(-1,0),(0,  t),  (1,0)}.  Then  H{0)  =  lim(_o  H{t)  and,  for  t  5^  0,  H{t)i  =  tti, 
while  ff(0)i  =  ir2(lR)  o  ()'’°.  We  will  show  that  this  example  of  a  discontinuity  is  prototypical. 

For  want  of  a  better  word,  we  call  the  linear  subspace  H  of  Ao  regular  in  case  is  of 
least  degree,  i.e.,  dim  Hi  n  7r,„  =  max{dim  H"! ,  dim  7r„i }  for  all  m.  Equivalently,  H  is  regular  iff 
HiHiTm  =  Hi  or  Xm,  for  all  m.  Thus,  H  is  regular  iff,  for  some  m,  7r<,„  Q  Hi  Ci^m-  We  note  in 
passing  that  this  m  provides  the  local  approximation  order  of  H . 

The  importance  of  this  notion  of  regularity  for  the  continuity  of  the  map  H  ^  Hi  is  illustrated 
by  the  following. 
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(6.2)  Lemma.  If  each  H{t)  is  re'  iJar  and  P  =  exists,  then  P  is  regular. 

Proof:  Since  H{t)i  converges,  its  dimension  is  eventually  constant,  and,  since  each  is  regular,  this 
implies  the  existence  of  some  m  so  that,  for  all  small  t,  each  H{t)i  is  scale-invariant  and  satisfies 
^<m  C  n{t)i  C  TTm,  hence  their  limit  P  satisfies  the  same  conditions.  4 

(6.3)  Lemma.  The  Unite- dimensional  subspace  H  is  regular  if  and  only  if,  for  some  regular, 
scale-invariant  polynomial  space  P,  P  is  minimally  total  for  H. 

Proof:  Since  =  Hi,  H  is  regular  iff  Hi  is  regular,  hence  (4.3)Theorem  proves  the  “only 

if”. 

For  the  converse,  let  P  be  a  regular  scale-invariant  polynomial  space.  Then,  with  m  such  that 
dimjr<m  <  dimP  <  dirnTr^,  we  must  have  T<m  C  P^  =  P  C  Assume  that  P  is  minimally 
total  for  H.  Then  H  contains  a  dual  set  for  with  Pa  ■=  ()".  Consequently, 

for  all  )a|,l/3|  <  m,  D^ha{Q)  =  CaSap  for  some  Ca  ^  0,  hence  hai  =  Cc»()“-  This  implies  that 
7r<m  C  Hi-  To  see  that  Hi  C  Xm,  observe  that,  since  P  is  total  for  H,  we  can  find,  for  any 
h  €  P^\0,  some  p  G  P  so  that  p’h  0,  and,  since  P  C  Xm,  this  implies  that  hj  G  7r„,.  ^ 

(6.4)  Theorem.  The  set  of  regular  (finite-dimensional)  subspaces  of  Aq  is  open  and  dense  (in  the 
space  of  all  finite-dimensional  subspaces,  and  in  the  “gap  topology”). 

Proof:  Assume  that  H  is  regular.  We  have  to  prove  that  all  K  in  some  neighborhood  of  H  are 
also  regular.  But  this  follows  from  (6.3)Lemma  and  from  the  fact  that  if  A  is  minimally  total  for 
H ,  then  it  is  minimally  total  for  all  nearby  K.  This  proves  that  the  set  of  regular  H  is  open. 

To  prove  that  the  set  of  regular  subspaces  is  dense,  let  H  be  an  arbitrary  finite-dimensional 
subspace  and  let  P  be  an  arbitrary  regular  scale-invariant  polynomial  space  of  the  same  dimension 
as  H.  Consider  the  map  R  :  P  -*  H'  :  p  >-*  If  P  is  1-1,  H  itself  is  regular,  by  (6.3)Lemma. 
Otherwise,  choose  an  orthonormal  basis  (rj)"  for  P  whose  first  m  terms  span  kerP.  Then  H 
contains  a  dual  set  for  and  this  can  be  extended  to  a  basis  (hj)”  for  H.  Define 

kj  :=  hj  -f-  crj,  all  j.  Then 

{cSij,  i,j  <  m; 

0,  i<m<j; 

(l-\-€)6ij,  i,j>m, 

i.e.,  the  matrix  is  lower  block  triangular  with  diagonal  blocks  dm  and  (1  -|-  €)Im-n, 

hence  invertible  for  all  positive  e.  Consequently,  P  is  minimally  total  for  K  span  (hj  A  erj)", 
for  all  positive  e.  This  shows,  with  (6.3)Lemma,  ihat  every  neighborhood  of  H  contains  regular 
subspaces.  ♦ 

Remark.  The  proof  actually  shows  that,  in  the  set  of  all  rz-dimensional  subspaces  of  Aq,  those  for 
which  a  fixed  regular  polynomial  space  is  minimally  total  form  an  open  and  dense  set. 

(6.5)  Corollary.  The  map  H  >->■  Hi  is  continuous  at  H  if  and  only  if  H  is  regular. 

Proof;  Assume  that  H  Hi  is  continuous  at  P.  By  (6.4)Theorem,  we  can  find  regular  H(t) 
with  limt_o  H(t)  =  H,  hence,  by  the  a.ssumed  continuity.  Hi  =  lim(_o  R(t)i,  thus  H  is  regular  by 
(6.2)Lemma. 

For  the  proof  of  the  converse,  assume,  more  precisely,  that  7r<fc  C  C  tt*.  By  (3.2)Lemma, 
this  implies  the  e.xistence  of  a  continuous  linear  projector  Tk,H  on  Ao  onto  P  which  preserves 
Taylor  polynomials  of  degree  <  k,  i.e.,  which  satisfies  TkTk,H  =  Tk-  Since  Tk,H  is  a  continuous 
linear  projector,  it  carries  any  subspace  K  sufficiently  close  to  P  1-1  onto  P.  Consequently,  for 
every  f  6  H,  we  can  find  in  every  such  K  an  element  g  close  to  /  which  satisfies  Tkg  =  Tt/.  In 
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particular,  =  /|  in  case  degfi  <  k,  hence  7r<fc  C  Ki  for  any  such  K.  Further,  if  deg/^  =  fc,  then 
Tkg  =  0,  i.e.,  deg^^  >  k,  hence,  since  g  is  close  to  /,  degfif^  =  k  and  gi  is  close  to  /|.  This  shows 
that,  for  such  K  and  every  p  6  Hi,  there  is  some  q  G  Ki  close  to  it,  and,  since  dim/ifi  =  dim  Hi 
for  such  K,  it  follows  that  Ki  and  Hi  are  also  close.  4 

7.  Hermite-BirkhofF  interpolation 

In  this  section,  we  consider  interpolation  by  polynomials  using  interpolation  conditions  of  the 
form  \0\p{D),  with  [B\  the  linear  functional  of  point  evaluation  at  6  and  p  a  polynomial.  More 
precisely,  we  want  to  interpolate  from  some  polynomial  space  Q,  using  the  interpolation  conditions 

(7.1)  A  :=  A(0;(Ptf))  :=  Y}^]Pe{D)  :=  :p  e  P,}. 

$  $ 


For  the  analysis  of  this  problem,  observe  that,  in  terms  of  (4.2), 

(7.2)  [^\p{D)q  =  p'E^q  =  q‘’(€sp),  Vp,  q  E 

This  implies  that  our  interpolation  problem,  as  specified  by  Q  and  A  =  A(Q;(P$)),  is  correct  if 
and  only  if  the  dual  problem  of  interpolation  from  H  :=  ^ePe  with  interpolation  conditions 

Q*  is  correct.  Therefore,  (4.5)Theorem  provides  the  following  conclusions. 

(7.3)  Theorem.  Given  any  finite  set  0  C  C*  and  corresponding  finite-dimensional  polyno¬ 
mial  spaces  Pg  {or  6  £  Q,  let  H  :=  '^gesPe-  Then  Hi  is  a  polynomial  space  of  least  degree 
among  all  polynomial  spaces  from  which  interpolation  with  interpolation  conditions  A(0;  (Ps))  := 
span{[0]p(i))  :  p  E  Pg;6  £  0}  is  correct. 

In  univariate  Hermite-Birkhoff  interpolation  [LJR],  one  matches  certain  derivatives  rather  than 
linear  combinations  of  derivatives.  Correspondingly,  we  will  use  the  term  Hermite-Birkhoff  in¬ 
terpolation  in  the  multivariate  context  of  (7.3)Theorem  in  case  all  the  spaces  Pa  have  a  homo¬ 
geneous  basis,  i.e.,  are  scale-invariant.  If  all  the  spaces  Pg  are,  in  addition,  /^-invariant,  then  we 
speak  of  Hermite  interpolation. 

If  each  Pg  is  P-invariant,  then  so  is  H  ;=  ^gCgPa,  therefore,  by  (3.12),  so  is  Hi-  In  the 
univariate  case,  it  follows  that  Hi  =  iTfc  for  some  k,  since  tt*  is  the  only  P-invariant  polynomial 
space  of  dimension  fc  -f  1.  Thus  we  obtain  the  wellknown  fact  that  univariate  Hermite  interpolation 
from  TTfc  with  fc  -f  1  conditions  is  always  correct. 

8.  Lagrange  interpolation 


The  special  case 

0  C  IR^  Pa  =  TTcVfi 

in  (7.3)Theorem  is  particularly  striking.  The  claim  here  is  that,  with 

expQ  :=  span(e9)se©, 

the  polynomial  space 

TT©  :=  (expe); 

is  of  least  degree  among  all  those  from  which  interpolation  at  0  is  uniquely  possible.  (We  are  using 
the  fact  that,  for  this  case,  H  =  H.)  See  (3.14)  for  a  recipe  for  generating  jt©. 


(8.1)  Example.  As  a  simple  illustration,  consider  s  =  2. 

For  #0  =  1,  TTe  =  xq. 

For  #0  =  2,  X0  =  7ri(IR)  o  (A-),  with  A  any  nonzero  vector  parcdlel  to  the  affine  hull  of  0. 
For  #0  =  3,  TTe  =  7r2(lR)  o  (A-),  with  A  any  nonzero  vector  parallel  to  the  affine  hull  of  0,  in 
case  that  hull  is  a  line.  Otherwise,  ire  =  xi. 

For  ^0  =  4,  xe  =  7r3(]R)  o  (A-),  with  A  any  nonzero  vector  parallel  to  the  affine  hull  of  0,  in 
case  that  hull  is  a  line.  Otherwise,  Xi  C  TTe  C  7r2-  In  that  case,  we  can  compute  the  barycentric 
coordinates  of  one  point  with  respect  to  the  other  three,  say  64  =  with  “(j)  =  I- 

On  setting  o(4)  =  -1,  we  thereby  obtain  the  (essentially  unique)  quadratic  homogeneous  element 
Note  that  its  span  is  a  continuous  function  of  0  except  when  the  four  points 
become  coUinear.  since  exp©  fails  to  be  regular  only  in  that  case. 

For  the  particular  choice  0  =  {O,0,r,0  +  r}  with  0  =  ei,r  =  ae2,  the  quadratic  term  becomes 

((er)^  +  (oe2-)^  -  ((^i  + 

hence  ttq  =  xi^i,  the  space  of  bilinear  polynomials.  Correspondingly,  the  interpolation  at  0  from 
X0  is,  in  this  case,  the  tensor  product  of  linear  interpolation,  as  one  would  hope.  4 

We  now  consider  how  x©  changes  with  0.  Since  e^+a  =  we  conclude  from  (3.7)Proposi- 
tion  that 

(8.2)  ^©+a  =  T^Q,  Va  6  IR. 

Further,  since  the  limit  at  0  of  any  H  C  Aq  does  not  change  under  scaling, 

(8.3)  iTa©  =  TT©,  Va  €  IR\0. 

Both  of  these  facts  could  also  be  deduced  directly  from  (3.14).  More  generally,  if  A  is  any  matrix, 
then  e^e  =  eg  0  A^,  hence,  by  (2.3), 

Xy4©  ^0  ®  A 

for  any  invertible  matrix  A. 

More  substantial  changes  in  0  may  change  x©  substantially.  In  fact,  the  map  0  i-+  x©  has 
jumps,  as  can  be  expected  from  (6. 5) Corollary.  The  simplest  possible  example  occurs  with  s  =  2 
and  #0  =  3  (cf.  (8.1)).  Here  x©  =  xi  except  when  0  is  collinear,  in  which  case  x©  =  xi(lR)  o  (A-) 
for  any  A  not  zero  and  parallel  to  the  affine  hull  of  0.  In  terms  of  the  Gram-Schmidt  algorithm 

(5.4)  foi  the  construction  of  x©  =  (exp0)|,  these  jumps  are  due  to  the  fact  that,  as  0  changes,  the 
degree  of  some  qk^  may  jump  even  if  we  arrange  the  starting  basis  (pj)  for  x©  in  such  a  fashion 
that  the  degree  of  each  q*,  ^  is  as  small  as  possible. 

Connected  with  this  is  the  fact  that,  near  a  0  at  which  the  map  jumps,  our  interpolation 
scheme  is  badly  behaved.  Put  positively,  in  that  caise,  it  is  a  much  stabler  thing  to  use  x©  when 
matching  data  at  some  0'  near  such  0,  rather  than  using  x©<.  Thus  in  our  simple  example,  it  would 
be  better  to  interpolate  at  0  from  Xi(IR)o(A-)  in  case  the  points  in  0  are  ‘nearly  collinear’,  rather 
than  from  x©  itself.  In  terms  of  the  Gram-Schmidt  algorithm  (5.4),  one  would  reject  a  ‘near-zero’ 
least  term  in  favor  of  the  next  non-zero  homogeneous  term. 

It  seems  natural  to  use  (e^jeg©  in  the  role  of  the  initial  basis  (pj)  in  the  Gram-Schmidt 
algorithm  (5.4),  making  use  of  the  fact  that  only  the  first  few  terms  in  the  Taylor  expansion  of  pk 
are  needed. 

We  can  also  use  induction  to  construct  the  unique  lef  ^  which  agrees  with  /  at  0,  thereby 
obtaining  its  Newton  form.  For  r  ^  0,  this  gives 

r  r  r  r  I  n  r  v 

f©Ur/  -  /©/  +  7- - nr7Z\(I  “ 

(1 
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with  Tr  :=  (ct  —  TeSr)  ^  the  least  term  of  the  error  when  interpolating  Cr  from  exp0  with  respect 
to  (ffe)*.  For,  Tt  G  jrauTN^re  by  (5.1)Proposition. 


9.  Osculation  and  coalescence 

Our  interpolation  scheme  Iq  depends  on  0  nicely  enough  to  allow  for  the  existence  of  a  limit 
when  some  or  all  of  the  points  in  0  coalesce  in  a  nice  enough  manner.  In  that  case,  the  limiting 
situation  often  is  Hermite  interpolation,  in  the  sense  defined  in  Section  7. 

Given  that  7©  is  characterized  by  t©  and  exp©,  it  is  natural,  in  light  of  the  remarks  following 
(4.1)Lemma,  to  study  the  limiting  situation  by  considering  the  limits  (if  any)  of  these  two  spaces, 
as  0  approaches  some  limiting  set  T. 

In  the  univariate  case,  ir©  =  7r<#©,  hence  it  does  not  change  (assuming  that  0  converges). 
Further,  exp©  converges  to 

®xpt,#t  •=  gT^<#-f 

t€T 

with  the  multiplicity  of  r,  i.e.,  the  number  of  points  from  0  which  coalesce  at  r.  In  particular, 
the  limit  always  exists  and  does  not  at  all  depend  on  just  how  0  approaches  T. 

The  multivariate  situation  is  much  more  complicated.  Neither  x©  nor  exp©  need  to  converge. 
If,  for  example,  s  =  2  and  0  =  {0,^}  and  6  alternates  between  the  two  axes  as  it  approaches  0,  then 
X©  alternates  between  the  span  of  ()°,  ()^’°  and  the  span  of  ()®,  ()°’^,  hence  does  not  converge.  Even 
if  X©  and  exp©  converge,  their  limits  strongly  depend  on  the  manner  in  which  0  approaches  T.  If, 
for  example,  6  =  6' t  +  o(t)  in  the  earlier  example,  then  limt_o  x©  =  Xi(I[l)  o  {9'-)  =  limt_o  exp©. 

Here  is  a  more  striking  example.  We  take  again  s  =  2,  but  talce  0  =  0(t)  =  {0,(t,0),(t^,t^)}. 
Then  0  is  in  general  position,  hence  x©  =  xi,  therefore  also  limt_o  5i‘©(j)  =  xj.  Further,  each  exp© 
contains  ()°.  In  addition,  exp©  contains  t"^(e(t.o)  ~  1)  *7Zo^  ()^’°*  Finely,  exp©  contains 

-  te(t,o)  -  (1  -  t))  ^  Q  '  0°’^  - 


Since  exp©  is  three-dimensional,  this  shows  that 

l^exp©  =  span  (()°,()^’°,()°’^  -  ()^’°)  ^  =  limx©. 

In  particular,  limexp©  is  not  even  scale-invariant,  hence  the  resulting  limiting  interpolation  scheme 
is  not  Hermite  or  Hermite-Birkhoff  interpolation  by  our  earlier  definition.  Also,  Urn  x©  fails  to 
equal  (limexp©)^.  But  limexp©  is  Z?-invariant  (as  it  has  to  be  since  each  exp©  is  7?-invariant). 

We  postpone  a  full  discussion  of  the  general  situation  to  a  future  paper  and  content  ourselves 
here  with  the  following  very  simple  case. 

We  assume  that,  more  precisely,  0  =  0(t)  consists  of  the  points  9{t),  with  9{t)  =  0(0)  9'(0)t. 

Consider  first  the  special  situation  that  0(0)  consists  of  one  point  only.  Then,  in  considering 
limt_o  exp©(t),  we  may  as  well  assume  that  0(0)  =  {0}  (see  (8.2)),  hence  0(t)  =  Et,  with  E  := 
{9'(0)  :  9  G  0}.  Consequently,  from  (6.1),  x=  =  (exps)^  =  limt_oexp=t,  while  This 

implies 

(9.1)  limx©(f)  =  XH  =  limexp©(t) 

in  this  simple  situation. 

The  same  argument  handles  the  slightly  more  general  situation  described  in  the  following 
proposition. 
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(9.2)  Proposition.  Assume  that  each  ff  E  O  =  Q(t)  is  of  the  form  Oq  +  0it.  Then,  for  each 
r  E  T  :=  {0o  ■  ^  G  0},  the  set  Hr  :=  {^i  :  Oq  =  t}  has  as  many  elements  as  there  are  0  G  0  with 
ffo  =  T,  and  limexpe  =  ICtct 

Proof:  The  exponential  space  exp©  is  the  direct  sum  of  the  spaces  exp^^t  and  these  converge 
to  er^TE,,  by  the  earlier  argument.  4 

(9.3)  Corollary.  If,  in  addition,  limexpe  is  regular,  then  Iq  converges  (e.g.,  pointwise  on  w)  to 
Hermite  interpolation  from  P  :=  (limexpe)^  with  interpolation  conditions 

Proof:  From  (6.5)Corollary,  we  conclude  that  the  assumed  regularity  of  lim  exp@  ensures  the 
convergence  of  ttq  to  P,  and  the  rest  follows  from  the  remarks  following  (4.1)Lemma.  4 

We  have  concluded  from  (6.1)  that  the  exponential  space  exp^©  approaches  the  polynomial 
space  X©  as  t  -+  0.  This  polynomial  space  is  finite-dimensional,  scale-invariant,  and  i?-invariant. 
It  would  be  very  nice  to  know  whether  every  finite-dimensional  scale-  and  D-invariant  polynomial 
space  arises  in  this  way.  For  it  would  then  be  possible  to  view  all  (regular)  osculation  as  coalescence. 
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